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We develop a Landau Ginzburg theory of the hidden order phase and the local moment antifer- 
romagnetic phase of URu2Si2. We unify the two broken symmetries in a common complex order 
parameter and derive many experimentally relevant consequences such as the topology of the phase 
diagram in magnetic field and pressure. The theory accounts for the appearance of a moment under 
application of stress and the thermal expansion anomaly across the phase transitions. It identifies 
the low energy mode which is seen in the hidden order phase near the conmensurate wavector (0,0, 
1) as the pseudo-Goldstone mode of the approximate U(l) symmetry. 



URii2Si2 is arguably the most intriguing heavy fermion 
material and its electronic structure has continued to be 
the focus of intensive investigations. It displays a phase 
transition at low temperatures material to a phase for 
which, in spite of a large number of experimental ef- 
forts, the order parameter has not been identified and 
it is therefore referred to as hidden order (HO) pQ. 

Recent first principles LDA+DMFT calculations 
showed that in the paramagnetic phase of URu2Si2, the 
local ground state of the configuration, and its first 
excited state are two singlets separated by a crystal field 
splitting A. In an atomic picture, the corresponding low 
lying singlets are: 

10} = 7l(l 4 ) - I - 4)), and |1> = ^(| 4 ) + | - 4)) - 
sin(0)|O), with <j> ~ 0.372tt and A w 35K. 

We proposed that the order parameter of the hidden 
order (HO) and local moment antiferromagnetic (LM A) 
phase is the excitonic mixing between the two lowest ly- 
ing configurations of the f electrons [2 . The excitonic 
mixing of the two singlets is described by a complex or- 
der parameter (X$i(j)) = ipj/2 = (ipij +iip2,j)/2, where 
-^01 (i) is the Hubbard operator |0)(1| at site j. ip2 is 
proportional to the magnetization along the z axis in the 
material, while ip\ is proportional to the hexadecapole 
operator (J x J y -\-J y J x )(J^ — Jy). The characteristic shape 
of the hexadecapole is shown in Fig. [TJ^. The presence 
of ifji does not break the time reversal symmetry, nor 
the tetragonal symmetry, but it does break the reflec- 
tion symmetry along the x and y axis. We identify the 
phase with nonzero ipi as the "hidden order" phase, and 
the phase with nonzero ^2 as the LMA phase. In this 
picture, the LMA and the HO order parameters are in- 
timately connected, and they are related by an internal 
rotation in parameter space. 

In this publication, we build on these insights from 
LDA+DMFT microscopic theory to construct a low en- 
ergy model, and to establish contact with many of the 
available experimental results on this material. Phe- 
nomenological theories of the Landau Ginzburg type for 
URu2Si2 have been developed before [HI SHU- However 
an approximate symmetry between the LMA and and 
the HO phase has not been noticed before, and the ma- 



terial specific information resulting from the microscopic 
calculations, was not available before. The new insights 
from LDA+DMFT calculation restricts the effective the- 
ory enough to result in a large number of consequences 
that can be compared with experiment, as for example 
the response of the system to pressure, uniaxial stresses 
and weak external magnetic field. In addition, the simpli- 
fications offered by the low energy effective Hamiltonian 
of Landau- Gindzburg type allow us to obtain analytical 
expressions which are not possible to obtain in the full 
LDA+DMFT solution. 

It is illuminating to interpret the two low lying con- 
figurations of the U atom as a local two level systems 
of (pseudo) spins 1/2 at each lattice site. The cor- 
respondence involves the three Pauli matrices: cr 3 = 
|0><0| - |1><1|, *i = |0}(1| + |1>(0|, a 2 = i(|l)(0| - |0)(1|). 
The crystal field splitting A of the two singlet states 
plays the role of external magnetic field in the third di- 
rection, while the real and imaginary part of the order 
parameter appear as the magnetization along the first 
and second direction in an internal order parameter space 
((01) = ^1? ( a 2) = ^2)- The coupling to the external 
magnetic field is given by B\±b{L z + 2S Z ) and since only 
the off-diagonal terms are nonzero (1|J^|0) = 4zcos((/>) 
and they are purely imaginary, the coupling is propor- 
tional to |(l|Z,s + 2S z |0)|£<7 2 . 

The coupling between uranium atoms is here modeled 
by a set of exchange constants Jg, which we allow to 
be different in the two different ordered states a = 1,2. 
Ignoring for the time being the coupling to the fermionic 
quasiparticle excitations, the low energy effective theory 
can be described by 

H = ^~*i-ii B B\(l\L z + 2S z m<r2 

i 

+ Y.V J 'A°{ + J iA4)- a) 

This Hamiltonian should be regarded as a simplified toy 
model to capture symmetry related changes not to far 
from the phase transitions. The LDA+DMFT treat- 
ment indicates that a more refined treatment should in- 
clude the coupling to dispersing electronic excitations. 
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This coupling is important for understanding several 
properties, including the low tempeature superconduc- 
tivity and the fermi surface reconstruction as a result 
of the coupling of the fermions to the excitonic mode. 
The LDA+DMFT Fermi surface has been shown to 
give rise to partially nested features with a wavevec- 
tor (0.6, 0, 0) [2]. Consequently, the itinerant electrons 
control the mangetic response around this wavevetor as 
shown in Ref. 6 . The analysis of the coupling of the mode 
to the fermionic quasiparticles is left for future study. 

We now treat this Hamiltonian in the mean field ap- 
proximation and arrive at an effective free energy of the 
spin system, written in terms of the order parameters 
ipaj and the conjugate Weiss fields h a ^ [7] 



where = EjVw) 2 and 5 = f coth(/?A/2) while 



[sinh(/?A) 



/ o A] cosh^A/2) 
' sinh 4 (/3A/2) * 



We determine the effective exchange constants J e f f at 
zero pressure in such a way to reproduce the experimen- 
tally observed critical temperatures = 

A 



and Jl ff - tanh ( A /(2T iV )) 



tanh(A/(2T )) 

where T = 17.7 K is the hid- 
den order transition and T/v = 15.7 K is the Neel tem- 
perature. The exchange constants are strain dependent 
and for the small compression we may use the expansion 



ja 



xx + e yy)) 



F[h, tj)] = - ^2 J ij ^i^aj ~ ( ftc M + b a )^a 

ij,a=(l,2) i,a=(l,2) 

-^T^log (cosh (py/(A/2) 2 + (h lti ) 2 + (h2 ti ) 



Here b 2 = b = Bfi B \(l\L z + 2^|0)|, with \{1\L Z + 
25^|0)| = 3.2cos((/>) « 1.25, is proportional to the exter- 
nal magnetic field B, while b\ is a fictitious field which 
couples to the hexadecapole order. The imaginary part 
of ip breaks the time reversal symmetry and therefore 
couples to the magnetic field. The mean field equations 
can be obtained by extremizing the free energy Eq. (|| 



h a j tanh(/3A^) 
"~2 



where e xx and e yy is the strain (compression) in x and y 
direction. The two constants g\ and g 2 are determined 
in such a way to reproduce the experimental variation 
of transition temperature with pressure T c (p). We take 
Si = 20, g 2 = 49. 

Having determined all the parameters of the Landau 
(2) Ginzburg functional, we proceed to determine the phase 
diagram as a function of pressure and magnetic field, 
which can be directly compared to the recent experimen- 
tal results of Refs. [ID|ITT]. 

There are two solutions of the mean field theory, the 
local moment antiferromagnetic solution with nonzero 
staggered ip2,i and vanishing ipi^. The second solution 
has nonzero staggered ipi^ but vansishing ip2,i, hence it 
shows no magnetic moment and corresponds to the hid- 
den order phase. 



(3) 
(4) 



with A* = ^/(a/2) z + (h hi y + (h 2ti y. 

We take the exchange constants between uranium sites 
in body centered tetragonal structure to be ferromagnetic 
in the same plane (RKKY is ferromagnetic at short dis- 
tance), and antiferromagnetic in c direction, favoring the 
staggered order with wave vector Q = (0,0, 1), observed 
experimentally [5J |9j [10] . As a result of the close simi- 
larity of the exchange constants of the hexadecapole and 
the LMA order parameters, both condense at the order- 
ing vector Q = (0, 0, 1). 

When the restriction to short range exchange is im- 
posed, the only combination of the exchange constants 
that enters the mean field equation is J e ff = 4| J\ \ +8J2, 
where J\ is the nearest neighbor ferromagnetic exchange 
and J2 is the antiferromagnetic exchange in c direction. 
The critical temperature in the absence of magnetic field 
is then given by T c = A/ (2 atanh(A/J e //)). 

Near the transition, when the field ip is small, the free 
energy acquires a simple Gindzburg-Landau form 




Hexadecapole 



FIG. 1: (a) The phase diagram of the mean field theory 
under pressure and in magnetic field (b) Experimental phase 
diagram by Aoki et al 11 . (c) The symmetry of the hex- 
adecapole order parameter of the HO phase. 



ij,a=(l,2) 



bifj 2 ,i Figure [T^i shows the phase diagram of the Ginzburg- 
Landau theory in the magnetic field under applied 
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pressure and Fig. [T}d show the experimentally de- 
termined phase diagram [11]. At low pressure, 
with decreasing temperature, there is a second or- 
der phase transition into the hidden order state 
at temperature Tq = wzzzzzzrS^ , T i — \ with \ c 



2arctanh(AA c / J\ ff ) 



1 + (#) 2 ( n Je i J 72 ) • At a critical field, b c = 

\ J eff + J e ff J 

f (^ff + J e 2 // )^/l/A 2 -l/(J e 1 // ) 2 , the hidden order 
phase is replaced by the fully polarized paramagnetic 
phase. 

At higher pressures, there is a second order tran- 
sition into the antiferromagnetic state at temperature 

A v /l+(26) 2 

1m — v -j — - — . 1 he two phases are sepa- 

2arctanh(Ayi+(26) 2 /j2 // ) * * 

rated by a first order boundary shown in Fig.[]Ji. It is also 
apparent from the figure that magnetic field under high 
pressure stabilize the hidden order phase relatively to the 
antiferromagnetic phase. The hidden order phase in a 
magnetic field develops a uniform magnetization directly 
proportional to the magnetic field M = Bfj^\(l\L z + 
2^|0)| 2 |/(J e 1 // + J 2 eff ) « 0.01/i B (£/T), in good quanti- 
tative agreement with experiment of Ref. [T2j Our model 
also allows us to extract the jump of the specific heat 
across the hidden order to paramagnetic phase bound- 

ary. It is given by Ac v = (a) 

hence Ac v /T c w 245 mJ/mol K 2 , and is in good agree- 
ment with experiment [T3] . 
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FIG. 2: The phase diagram (a) under pressure and (b),(d) 
stress in ab plane, and uniaxial strain (c). The critical tem- 
peratures below which the HO and the LMA solutions are 
possible are indicated in the diagram. The regions where 
LMA and HO are stable are shaded in blue and green respec- 
tively. The transition between the HO and the LMA phase is 
of the first order. 

The pressure dependence of the critical temperature 
is shown in Fig. [2^i. The close similarity with experi- 
ment of Hassinger et al. [14] is expected, since the the 
couplings gi and #2 were optimized to reproduce correct 
pressure dependence of T c . However, strain in a and 



c direction acts in quite unexpected way on the stabil- 
ity of the two phases. The stress is connected to strain 
by the elastic constants (a = Ce)(see Ref. [16] for de- 
tails). Their value was determined by the ultrasonic- 
sound-velocity measurements to be (en, C33, C12, C13) = 
(255,313,48,86) x 10 10 erg/cm 3 Q2|. 

Application of strain in ab plane (<r|| (1,0,0) or 
<j\ |(1, 1,0)) destabilizes the hidden order state, and sta- 
bilizes the antiferromagnetic state at a around 0.6 GPa. 
Within our model, the transition temperatures for the 
two phases are given by 



atanh 



A/2 



J eff^ 



9a (c 11 +C 12 C )233-2c2 3 Cr ^) / 



Hence LMA is nucleated upon application of very small 
in-plane stress. This has been observed in the neutron 
scattering experiments of Ref. 16. In those experiments, 
the staggered magnetization was found to grow linearly 
with the applied stress, an observation which has been 
taken as evidence for time reversal breaking in the hid- 
den order phase [18]. In our framework, these measure- 
ments [16] can be interpreted as the result of inhomo- 
geneities in the strain field, which given the very low 
barrier between HO and the LMA state, easily nucleates 
LMA regions. The LMA phase persists as a metastable 
state all the way to zero strain, indicating that nucleation 
of LMA regions is actually possible for infinitesimal stress 
in the presence of defects. This picture can be tested by 
scanning tunnelling microscopy of stressed samples. 

So far, we have considered stress breaking the tetrag- 
onal symmetry. Uniaxial strain stabilizes hidden order 
(the T^eei smaller than T ) and decreases the hidden 
order transition temperature, as shown in Fig. The 
formula for the two transition temperatures is given by 



T?(a zz ) = 



atanh 



A/2 



9a 



2ci3 



(cil+Ci2)c 3 3-2cf 3 



,) 



The 5.6 GPa strain leads to complete elimination of the 
hidden order. However, the uniaxial strain also affects 
the crystal field splitting A. We have not attempted 
to model this dependence quantitatively in our theory, 
but we notice that an increase in A will result in a 
rapid decrease of the hidden order temperature (and the 
LMA critical temperature) resulting in a quantum crit- 
ical point at critical value of the crystal field splitting 
A c = J 1 e ff separating the HO phase from the paramag- 
netic phase at zero temperature. This critical point may 
already have been accessed via Rh doping in Ref. [20] This 
quantum critical point is expected to be in the same uni- 
versality class as that of Ising itinerant antiferromagnets 

ED- 

To study the dilatation of the sample in a and c di- 
rection due to the emergence of the HO and LMA states 
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FIG. 3: The strain resulting from the phase transition into 
HO and LAM state. The sample contracts in ab plane (posi- 
tive strain), and expands in c direction (negative strain). 

we add to the Landau free energy the standard elasticity 
terms and their coupling to the order parameter 

^eCe -ea+^ ^ Jg(l + g a (e xx + e yy )) VwVw- 

zj,a=(l,2) 

Here a and e are the stress and strain tensor and 
the tensor C describe the elastic moduli. Differentiat- 
ing this free energy we obtain the dilatation in a and 
c direction, e xx = c^L a /2 and e zz = —cisL a , where 

^ (Jcil+C 12 )c33-2cl 3 ) ' 

Fig. [3] shows the temperature dependence of the strain 
in the a and c direction at zero pressure, where the tran- 
sition into the hidden order state occurs. The sample 
contracts in the afr-plane (positive strain) and expands 
in c direction (negative strain). The expansion in c di- 
rection is smaller than expansion in ab direction. At a 
pressure of 2 GPa, where the LMA state is stable, the di- 
latation has the same trend but the magnitudes are con- 
siderably larger. Similar temperature variation of the di- 
latation was recently measured in Ref. 22, where consid- 
erably larger dilatation was found in transition to LMA 
phase than to the hidden order phase. This difference 
is connected to the larger slope of the LMA transition 
temperature compared to hidden order transition tem- 
perature. 

The Landau Ginzburg free energy explicitly exhibits 
the remarkable similarity between the hidden order phase 
and LMA phase. There is an approximate U(l) rota- 
tional symmetry in an internal parameter space. The 
main difference between the two phases lies in the dif- 
ferent values of the exchange constants, which is of the 
order of 6 %. Our theory therefore provides a microscopic 
basis for the remarkable similarity between the LMA and 
HO phase, dubbed adiabatic continuity [T5] . 

If the U(l) symmetry was exact, the spontaneous 
breaking of this symmetry would result in a Goldstone 
mode describing the transverse fluctuations of the order 
parameter, i.e., the potential would have the form of a 



Mexican hat with a flat bottom. Due to the small ex- 
plicit breaking of the U(l) symmetry, since the exchange 
constants and J 2 ^ are slightly different, this Gold- 
stone mode acquires a small mass, and we refer to it as 
a pseudo-Goldstone mode. In the hidden order phase 
7^ the pseudo-Goldsone mode can be identified 
with ip2,i (the transverse fluctuation in the Mexican hat). 
Hence the pseudo-Goldstone mode of the hidden order 
phase carries a magnetic moment and can be observed 
by neutron scattering. In the antiferromagnetic phase the 
pseudo-Goldstone mode can be identify with ipij, which 
carries hexadecapolar moment but no mangetic moment 
and is invisible to neutrons. As a result a low lying mode 
at (0, 0, 1) is only visible to neutrons in the hidden order 
phase. This provides a natural explanation of the mode 
observed in neutron scattering experiments, as measured 
by Broholm et al [8j E] and Villaume et all [10 . The 
mode was observed only in the hidden order phase, but 
not in the antiferromagnetic phase [10] . The energy scale 
of this mode is a measure of how different the hidden or- 
der phase is from the antiferromagnetic phase, and there- 
fore should decrease with pressure (since the difference in 
the exchange constants decreases with increasing pres- 
sure) and should increase with magnetic field (since the 
magnetic field destabilizes the antiferromagnetic phase 
further relative to the hidden order phase). 

In conclusion, we developed an effective theory of 
the paramagnetic to hidden order and local antiferro- 
magnetic transition. The theory is consistent with a 
large body of experimental data, is inspired by micro- 
scopic LDA+DMFT calculations, and puts URu2Si2 in a 
broader context of other f 2 systems [23j |24] . 
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